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Abstract 

We derive the expression for the anisotropic heterogeneous matrices of per- 
mittivity and permeability associated with two-dimensional polygonal and star 
shaped cloaks. We numerically show using finite elements that the forward scat- 
tering worsens when we increase the number of sides in the latter cloaks, whereas 
it improves for the former ones. This antagonistic behavior is discussed using a 
rigorous asymptotic approach. We use a symmetry group theoretical approach 
to derive the cloaks design. 

Ocis: (000.3860) Mathematical methods in physics; (260.2110) Electromagnetic the- 
ory; (160.3918) Metamaterials; (160.1190) Anisotropic optical materials 

1 Introduction 

Back in 2006, a team led by Pendry and Smith theorized and subsequently experimen- 
tally validated that a finite size object surrounded by a spherical coating consisting of 
a metamaterial might become invisible for electromagnetic waves [U [2] . Many authors 
have since then dedicated a fast growing amount of work to the invisibility cloak- 
ing problem. However, there are alternative approaches, including inverse problems's 
methods [3] or plasmonic properties of negatively refracting index materials pfj. It is 
also possible to cloak arbitrarily shaped regions, and for cylindrical domains of arbi- 
trary cross-section, one can generalize pQ to linear radial geometric transformations 


t> = Ri{9) + r{R 2 {9) - R 1 {9))/R 2 {9) , < r < R 2 {6) , 
$> = , < 9 < 2tt , (1) 

x' 3 = x 3 , x 3 e It , 

where r', 9' and x' 3 are "radially contracted cylindrical coordinates" and (xi, £ 2 , ^3) 
is the Cartesian basis. This transformation maps the arbitrary domain Dr^q) onto 
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the coating D R2 i re \ \ D Rl r r ^ with Ri{r,9) < R 2 (r,9). In other words, if a source 
located in H 2 \ D^rfi) radiates in a vacuum, the electromagnetic field cannot reach 
the domain D Rl ( r g\ and therefore this region is a shelter for any object. Moreover, the 
transformation maps the field within the domain Dr^^) onto itself by the identity 
transformation. In this paper, we would like to investigate designs of polygonal and 




{12/5} {12/4} {12/3} {12/2} {12/5} 



Figure 1: (a) Triangle of vertices (0, 0), 9t) and (b i+ i, 9 i+ i) used in the triangulation 
of a polygon (see (b)) or a polygram (see (c)) mapped under Eq. ([1]) onto the trapezoid 
of vertices (aj,^), (bi,9i), (b i+ i,9 i+ i) and (a i+ i, 9 i+ i). The angles 9{ and the segments 
dij are as in Eq. ([2]); (b) Dodecagon and first step in its four tessellations {12/2}- 
{12/5}; (c) Tessellation {12/5} and representative samples of its triangulation. 

star-shaped cloaks which we can triangulate, see Figs. [Q [2], [31 In this case, consider 
a given triangle, numbered i, of the triangulation of the region we want to cloak. We 
adapt Eqs. ([1]) to the specific geometry of triangles, to derive a geometric transform 
that maps the whole triangle of vertices (0,0), (^,6^), (6 i+1 ,6 l i+ i) onto a trapezoidal 
region of vertices (a.j,0j), (a i+1 ,9 i+ i), (6*4.1, 6> i+1 ), see Fig. [TJ This means that 

in Eqs. (CQ), the points (Ri(9),9) and (R2(9),9) which now respectively lie in the line 
segments [(a i; 9 { ), (a i+1 , 9 i+ i)} and 0*), (b i+1 , 9 i+ i)}, are referred to as (Ri tl (9), 9) and 
(Ri, 2 (9),9) where 

' Ri,j = dij cos(9 - 9 ), 0i < 9 <9 i+1 , j = 1,2 

dij = Cijc i+ld sin(6» i+ i - f )/ ^ c 2 ^ + c 2 +1J - 2c itj c i+ltj cos(9 i+1 - 9~), j = 1,2 
— 8i + &rccos(dij / Ci t j) = 9 i+ i + arccos(dj iJ /cj + ij) 

k c i,l = a i i c i,2 — 0, 
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In the sequel, we make use of symmetries to further simplify the implementation of 
such cloaks. More precisely, we work in a reference triangle and construct the overall 
cloak via compound reflection symmetries and rotations from this reference region. 



2 Change of coordinates and pullbacks from optical 
to physical space 

The basic principle of our approach is to transform a geometrical domain into another 
one and to search how the Maxwell equations to be solved have to be changed. As we 
start with a given set of equations on a given domain, it seems at first sight that we 
have to map this domain on a new one. Nevertheless, it is the opposite that has to 
be done. On the original domain, say Q, the equations are described in a particular 
coordinate system x, (taken here by default to be Cartesian coordinates). We want to 
establish a one-to-one correspondence with a new domain Q' where we use a general 
(possibly non-orthogonal) coordinate system yi. 

domain Q with < — domain Q' with 
coordinates Xj %i{Vj) coordinates yi 



differential dx« — > differential Y2j dy 
pullback 

More precisely, consider a map from a coordinate system {2/1,1/2,1/3} to the initial one 
{x 1 ,x 2 ,x 3 } given by the transformation x 1 (y 1 ,y 2 ,y 3 ), £2(2/1,2/2,2/3) and x 3 ( y y 1 ,y 2 ,y 3 )- 
Here, we must point out that, we are mapping the transformed domain and coordinate 
system onto the initial one with Cartesian coordinates. This change of coordinates is 
characterized by the transformation of the differentials through dx = 3 xy dy where 



d(xi,x 2 ,x 3 ) 
5(2/1,2/2,2/3) 



J xy = af .. ~ ? dx = ( da; i' dx 2, dx 3 ) T and dy = (dy 1 , dy 2 , dy 3 ) T . (3) 



In the sequel we will also need a compound transformation. Let us consider three co- 
ordinate systems {2/1,2/2,2/3}, { z i, z 2, z 3 }, and {xi,x 2 ,x 3 } (possibly on different regions 
of spaces) . The two successive changes of coordinates are given by the Jacobian ma- 
trices J xz and J zy . This rule naturally applies for an arbitrary number of coordinate 
systems and we call the resulting Jacobian matrix J = J xz 3 zy . In electromagnetism, 
changes of coordinates amount to replacing the different materials (often homogeneous 
and isotropic, which corresponds to the case of scalar piecewise constant permittivity 
e and permeability /i) by equivalent inhomogeneous anisotropic materials described by 
a transformation matrix T (metric tensor) [6J, [TJ, El E] so that the new permittivities 
and permeabilities read 

e!_ = eT- 1 , and // = //T^ 1 where T = J T J/det(J) . (4) 
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3 Transformation optics and the geometry of irreg- 
ular polygons 



3.1 Material properties of the regions 



From Eqs. (OQ), (J2J) and (HI), we derive the transformation matrix associated with a 
reference trapezoidal region (see Fig. [T|) used to generate, via reflections symmetry 
and rotations, two-dimensional polygonal and star-shaped cloaks: 
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(5) 



where 



en(6 , ) = R2(6')/(R 2 (6')-R 1 (6>)) 



(6) 



and 



(r> - R 2 (e'))R 2 (9')^P - (r> - R^R^e')^ 



(R 2 (6') - R 1 (9')y 



(7) 



for 



Ri(0') <r'<R 2 (e'), (8) 

with 

fr(r^) = (r'-R i m) R2 ^%, ) (9) 

(the expression for e\ 2 corrects that in [5]). Elsewhere, T _1 reduces to the identity 
matrix (en = 1, e\ 2 — and f r = r' for r' > R 2 (8')). However, Eqs. (jlj) and 
([5]) can just be applied to derive the (anisotropic heterogeneous) permittivity d_ and 
permeability // in the reference triangle. An elegant way to deduce optical properties 
of other regions of the triangulation is to use linear groups of symmetry in the plane, 
as now explained. 



3.2 Symmetries and transformation groups in cloaks' con- 
struction 

Here, we use Euclidean groups, namely reflections symmetry and rotations in the plane, 
as a tool to design very general polygonal and star-shaped cloaks. We choose our 
reference triangle (labeled as triangle number 1, or just referred to as first triangle) to 
be in the first quadrant with one of its sides belonging to the Xi-axis. In the regular 
polygons case, the triangle number % of the triangulation, is simply obtained from this 
first one, via the rotation with angle {% — 1)-. In the irregular case (star-shaped), due 
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to geometric chirality properties (irregular triangles are not identical to their mirror 
images), triangles of odd numbers cannot be obtained by rotations from the original 
one. More precisely, the neighboring triangles are obtained, from the first one, as mirror 
images with respect to the shared side (more precisely, a mirror image with respect to 
the plane perpendicular to the triangle and containing the shared side). Hence, any 
triangle of even number i, is obtained by applying to the first triangle the linear map 

71 71 71 71 

(xi, x 2 ) > (xi cos(— i) + x 2 sin(— i), x\ sin(— i) — x 2 cos(— i)) (10) 
n n n n 

which is the composition of the reflection (xx,^) i— > (x\, —x 2 ) with respect to xj-axis, 
followed by the rotation with angle -i. Whereas, every triangle of odd number i, is 
obtained by rotating the first triangle with angle — In short, optical properties of 
all regions within the triangulation are deduced from the reference trapezoidal region 
using compositions of rotations and reflections. This algorithm is particularly well- 
suited for all polygonal or star shaped regions displaying ra-fold symmetries. 

3.3 Tesselation and stellation 

We are interested in the antagonistic behaviour of polygonal and star shaped cloaks. 
There is a natural correspondence between these two types of cloaks, and this can be 
seen as follows. 

Consider a regular star polygon which is represented by its Schlafli symbol {n/m}, 
where n is the number of vertices, and m is the step used in sequencing the edges 
around it. For example, take a dodecagram {12/5} as in Fig. [T] which has twelve 
vertices, see Fig. [3]^c). One can draw a straight line joining each of its vertices to 
another one, which is five vertices apart. Or equivalently, choose a vertex, say, on the 
horizontal axis and project the other vertices onto that axis as depicted in Fig. [U(b). 
For a given dodecagon with Schlafli symbol {12}, there are four segments obtained by 
projection which generate four distinct stellations {12/2}-{12/5}. Note that a regular 
polygon with n sides has (n — 4)/2 stellations if n is even, and (n — 3)/2 stellations, 
if n is odd. We display in Fig. [2] and [3] a gallery of cloaks together with their Schlafli 
symbol. We note that the higher the Schlafli symbol of the polygons, the smoother the 
wave patterns (the less scattering). 

On the contrary, the scattering worsens for their corresponding stellations. This is 
not surprising since it is fairly intuitive that a cloak with many sharp angles will tend 
to display more singularities in the entries of the transformation matrix describing the 
electromagnetic material properties. However, we now give an asymptotic argument 
supporting this numerical observation. 

3.4 Asymptotic behaviour of the T -1 matrix 

The tessellation of regular polygons involves n regular triangles. We can chose 

71 71 

Oj = a i+ x =: Ri, h = b i+1 =: R 2 , 9 i+ i — 9i — — and 9 = s—, < s < 1. (11) 

n n 
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The coefficients of T now read 

IS 1 = (^V) 2r ^ + o(A), 

/t 2 — iti r rr 
rl _ ^ + 2^-^) 1 1 

Jl2 " " 2r'(r'-R 1 ) n +U( n3 j ' ^ 

and TjJ 1 = 0, otherwise. So as n grows, T -1 tends to 

= Dia ^^^' 7"^' (13) 



This is the expression of the transformation matrix for a circular cloak, see also [10] 
for a particular case. Now assume a i+1 ^ a i; or 7^ foj, which is expected to hold in 
star shaped case, the Taylor expansion shows that, at least one tends to infinity 
or det(T" 1 ) tends to zero, as n increases. For example, 

T ^ = TT1 V7 df±^^M±l n * + 0(1), (14) 

(r's{a i+1 - at) + a i+1 (r' - a i )){a i+1 + s{a i+1 - a*)) W 

if dj+i > a, and > &j. This explains the origin of the singularities in the material 
parameters of stellated cloaks. 



4 Finite elements computations 

Let us now compute the total electromagnetic field for a plane wave incident upon a 
cylindrical cloak. In p polarization, if e3 is a unit vector oriented along the axis of the 
cylinder, the longitudinal electric field E/ = E 3 (xi,X2)e 3 is solution of: 

V x (V _1 V x Ej) - k 2 £Ei = (15) 

where k = uj^/pqEq = u/c is the wavenumber, c being the speed of light in vacuum, 
and and p! are defined by Eqs. (BJ and (JSJ). Also, E; = Ej + E rf , where Ej = e lfcxi e 3 is 
the incident field and E^ is the diffracted field which satisfies the usual outgoing wave 
condition (to ensure existence and uniqueness of the solution). We have implemented 
the weak form of this scattering problem in the finite element package COMSOL (with 
Perfectly Matched Layers modelling the unbounded domain). We consider a plane wave 
incident from the left at wavelength A = 0.35 {all lengths are given in arbitrary units, 
pm for instance). We note that when we increase the number of sides of a polygonal 
cloak, the scattering becomes less apparent in Fig. [2], which is in agreement with the 
fact that the transformation matrix T _1 tends to that of a circular cloak, see formula 
(JTSJ). On the contrary, both forward and backward scattering worsen when we increase 
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the number of sides of corresponding stellated (star-shaped) cloaks, see Fig. [3J and we 
have identified the origin of this singular behaviour in the previous section. Note that 
we kept the frequency constant as well as the mesh (with around 60.000 elements), so 
as to exemplify that any attempt to manufacture a star-shaped cloak with many sides 
would be fairly challenging: when we refine the mesh, the scattering indeed reduces 
suggesting severer constraints on the corresponding metamaterial. Interestingly, the 
field can be confined when the frequency of the plane wave excites some resonance of 
the star-shaped cavity inside the cloak, see Fig. [3] (f). 

5 Conclusion 

Analysis of polygonal invisibility cloaks is of current interest [TUl [TT| [T2] . In this pa- 
per, we made use of a symmetry group action to construct such cloaks, and we have 
numerically shown that, for a given frequency, n-sided polygonal cloaks behave in a 
way similar to circular cloaks when n increases whereas their stellated counterparts 
become more and more singular. An asymptotic algorithm supports these observa- 
tions. A natural question to ask is whether it is possible to design such cloaks using a 
homogenization approach. Our analysis suggests that the answer should be positive for 
the former type of cloaks (as their boundary smoothens for large n and as it has been 
shown that structured circular cloaks can be homogenized [13]). We conjecture that 
the answer should be negative for the latter cloaks: it is well-known that homogeniza- 
tion theory fails in so-called bad domains that display boundaries with sharp corners 
(for which star-shaped invisibility cloaks is an electromagnetic paradigm). 
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